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Abstract 

We present a general procedure for determining possible (nonuniform) magnetic 
fields such that the Pauli equation becomes quasi-exactly solvable (QES) with an un- 
derlying sl{2) symmetry. This procedure makes full use of the close connection between 
QES systems and supersymmetry. Of the ten classes of ■s/(2)-based one-dimensional 
QES systems, we have found that nine classes allow such construction. 
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1 Introduction 



The Pauli equation describes the motion of a charged particle in an external magnetic field. 

It is given by the Hamiltonian {% = e = 2me = c = 1) 

H={p, + A,f + + Ayf + |(V X A),a, , (1) 

where p^;, Py are the momentum operators, 51 = 2 is the gyromagnetic ratio, A is the vector 
potential of the electromagnetic field, and cr^ is the Pauli matrix. For uniform magnetic field 
Bj. — By — 0, Bz — B the system is exactly solvable, giving the Landau levels. On the other 
hand, it was proved in [1] that for any general magnetic field B^ = B{x,y) perpendicular to 
the xy plane, the ground state is exactly calculable, owing to the existence of supersymmetry 
(SUSY) in eq.(l) [2]. The general result of [1] can be viewed as a very special case of the 
newly discovered quasi-exactly solvable (QES) systems, which are systems that allow parts 
of their spectrum to be solved algebraically ([4]-[13]). 

The Landau problem and the result of [1] represent the two extremes of the spectral 
problem of eq.(l). It is thus of interest to determine if other possibihties exist. Based on 
SUSY of the Pauli equation and the idea of shape invariance, it was shown that there exist 
three other forms of (nonuniform) magnetic field which make Pauli equation exactly solvable 
[2] (see also [3]). It seems difficult, if not impossible, to find other forms of the magnetic field 
such that the Pauli equation could be exactly solved. A more modest aim is to determine 
magnetic fields such that parts of the spectrum of the Pauli equation can be algebraically 
obtained. In other words, one looks for those magnetic fields under which the Pauli equation 
becomes quasi-exactly solvable (we note that for certain non uniform magnetic fields, the 
Pauh equation admits two solutions [14] and thus can be considered as a QES system). 
But even with this modest aim, the possibility seems enormous, since there are many QES 
systems based on different Lie algebras. In this paper we would like to make an attempt 
in this direction based on the simplest Lie algebra, namely, the sl{2) algebra. QES systems 
based on sl{2) algebra have been completely classified by Turbiner [5], and the necessary 
and sufficient conditions for the normalizability of the wave functions in such systems were 
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completely determined in [6]. It turns out that general forms of the magnetic field can be 
found so that the Pauli equation can be fitted into nine of the ten classes in [5] . The magnetic 
fields giving rise to these nine classes of QES potentials are divided into two groups: six in 
asymmetric gauge, and the other three in symmetric gauge. We shall describe these cases 
separately in Sect. 2 and 3. Wc would like to mention here that it is not necessary to 
consider magnetic fields that give rise to QES Pauli Hamiltonians with periodic potentials, 
since it has been proved in [6] that the wave functions in such systems are not normalizable. 



2 Magnetic field in asymmetric gauge 

Consider magnetic field in the asymmetric gauge given by the vector potential 

A,{x,y)^0 ,Ay{x,y)^-W{x) , (2) 

where W{x) is an arbitrary function of x. The magnetic field 

B has components — By — and B^ — B{x) = —W'{x). The Pauli Hamiltonian is 
then given by 

H = pl + {Py- Wix) f - W'ix)a, , (3) 

where W'{x) = dW{x)/dx. The eigenfunction ip can be factorized as 

iP{x,y) = e-'''yiP{x) . (4) 

Here k {—oo < k < oo) are the eigenvalues of Py, and i/jIx) is a two-component function of 
X. The upper and lower components of ip are then governed by the Hamiltonians and 
respectively, where 

= + {W{x) + k)\ W'{x) . (5) 

In this form the SUSY structure of the Pauli equation is clearly exhibited, with W{x) — 
W(x) + k playing the role of the supcrpotential. Once the upper component i/j oi i/j is 
solved for a nonzero energy, the lower component can be obtained by applying appropriate 
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supercharge on the upper component [2] and vice versa. Hence, it is suffice to consider the 
solution of the upper component, which satisfies the Schrodinger equation (we assume that 
the upper component has a normalizable zero energy state) H^%1) — Eip, where 

with 

V{x) = W{xy - W'{x) . (7) 

From the knowledge of shape invariant SUSY potentials, it was found that there are four 
allowed forms of shape invariant W{x) for which the spectrum of the Pauli equation can be 
algebraically written down [2] . One of the four forms gives rise to uniform magnetic field. 

The SUSY structure of the Pauli equation can be made use of in a different way, namely, 
in its close connection with quasi-exactly solvability [11, 12]. We shall give a brief description 
of this connection below. 

Consider a system described by eq.(6). We shall look for V{x) such that the system is 
QES. According to the theory of QES models, one first makes an imaginary gauge transfor- 
mation on the wave function ip{x) 

i;{x) = 0(x)e-^(^) , (8) 
where g{x) is the gauge function. The function satisfies 



dx'^ dx 

For physical systems which we are interested in, the phase factor ex^{—g{x)) is responsible 
for the asymptotic behaviors of the wave function so as to ensure normalizability. The 
function (p{x) satisfies a Schrodinger equation with a gauge transformed Hamiltonian 

where Wq{x) = g'{x). Now if V{x) is such that the quantal system is QES, that means 
the gauge transformed Hamiltonian Hq can be written as a quadratic combination of the 



generators J" of some Lie algebra with a finite dimensional representation. Within this finite 
dimensional Hilbcrt space the Hamiltonian Hq can be diagonalized, and therefore a finite 
number of eigenstates are solvable. For one-dimensional QES systems the most general Lie 
algebra is sl{2) ([4]- [11]). Hence if eq.(lO) is QES then it can be expressed as 

Hg^Y. CabJ^J^ + + constant , (11) 

where Cab, Ca are constant coefficients, and the J" are the generators of the Lie algebra 
sl{2) given by 

J+ = z'^-Nz, 
az 

J' = .^-f iV^0,l,2... 

= (^^) 

Here the variables x and z are related by 2; = h{x), where h{-) is some (explicit or implicit) 
function . The value j = N/2 is called the weight of the differential representation of sl{2) 
algebra, and N is the degree of the eigenf unctions, which are polynomials in a (A^ + 1)- 
dimensional Hilbert space. 

The requirement in eq.(ll) fixes V{x) and Wo{x), and Hq will have an algebraic sector 
with + 1 eigenvalues and eigenfunctions. In this sector the eigenfunction has the general 
form 

ijj ^ (^z - zi){z - Z2) ■ ■ -{z - ZN)ex.-p{- J Wo{x)dx) , (13) 

where Zi {i = 1,2, . . . , N) are N parameters that can be determined from the eigenvalue 
equations, namely, the Bethe ansatz equations corresponding to the QES problem [9, 13]. 
One can rewrite eq.(13) as 

iP = exp(- y" Wn{x, {z,})dx) , (14) 

and 



There are + 1 possible functions Wn{x, {zi}) for the + 1 sets of eigenfunctions tp. It is 

easy to check that satisfies the Ricatti equation [11, 12] 

- W = V - E , (16) 

where E is the cigenencrgy corresponding to the eigenfunction given in cq.(13) for a 
particular set of parameters {zi}. Eq.(16) shows the connection between SUSY and QES 
problems. 

Prom eqs.(6), (7) and (16) it is clear how one should proceed to determine the magnetic 
fields so that the Pauli equation becomes QES based on s/(2): one needs only to determine 
the superpotentials W{x) according to eq.(16) from the QES potentials V{x) classified by 
Turbiner [5]. This is easily done by observing that the supcrpotential Wq corresponding 
to = is related to the gauge function g{x) associated with a particular class of QES 
potential V{x) by g'{x) — Wo{x). Once Wo is obtained, then Bq — —Wq{x) is the required 
magnetic field that allows the weight zero {j = N = 0) state to be known in that class. But 
this state is just the ground state, and hence we have not gone beyond the result of [1] . 
What is more interesting is to obtain higher weight states (i.e. j > 0), which will include 
excited states. Por weight j {N = 2j) states, this is achieved by forming the supcrpotential 
Wn{x, {zi}) according to eq.(15). Of the 1 possible sets of solutions of the Bethe ansatz 
equations, the set of roots {zi, Z2, ■ ■ ■ , z^} to be used in eq.(15) is chosen to be the set for 
which the energy of the corresponding state is the lowest (usually it is the ground state). 
The required magnetic field which gives rise to the A'" + 1 solvable states is then obtained as 
Bn — —Wj^. Prom the table in [5] it is easily seen that only six classes need be considered, 
namely class I to class VI. Class VII to IX are excluded because these are QES systems with 
basic variables defined only on the half-line (0, oo), while class X corresponds to periodic 
potentials giving rise to non-normalizable wave functions. Below we shall illustrate our 
construction of QES magnetic fields through the class I and II QES systems, which serve as 
representative examples of two different types of QES problems. 
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Class I 

According to Turbiner's classification, the QES potential belonging to class 1 has the 
form [15] 

Vn{x) = a'e"'"^ - a [a{2N + 1) + 2bk] e""^ + c (26^ - a) e"" + c'e'"" + bl - 2ac . (17) 

Here bk = b + k with constant b. Without loss of generality, we assume a,a,c > for 
definiteness. The corresponding gauge function g{x) is given by 

g(x) = -e-"^ + -e"^ + bkX . (18) 
a a 

One should always keep in mind that the parameters selected must ensure convergence of 
the function exp{—g{x)) in order to guarantee normalizability of the wave function (this is 
generally not required by the mathematicians). We have also added the constant {bl — 2ac) 
in Vn so that for j — 0, the energy of the ground state is zero {E — Qi). This is not necessary, 
but it allows the results for j = and j > be presented in a unified way. The potential 
V[x) that gives the ground state is generated by 

V{x) ^ Vo-E 

= W^-W^, (19) 



with 



Wo{x) = g'ix) 

= -ae""^ + ce"^ + b + k . (20) 

The corresponding magnetic field is given by 

Bo = -W^^(x) 

= -aae-"^" - cae"^ . (21) 

To obtain magnetic fields and the corresponding potentials which admit solvable states 
with higher weights j, we must first derive the Bethe ansatz equations. To this end, let us 
perform the change of variable z = h{x) = exp{—ax). Eq.(9) then becomes 



2az' - {2bk + a)z - 2c 



dz 



+ 



E 

-2aNz 

a. 



, (22) 



which can be written as a quadratic combination of the sl{2) generators J"*", J and J° as 
^70 = 0, 

Ti = -Q;J+J- + 2aJ+-[a(iV+l) + 26fc] J°-2cJ" + constant . (23) 

For N > 0, there are + 1 solutions which include excited states. Assuming (f){z) = 
YliLii^ — Zi) in eq.(22), one obtains the Bethe ansatz equations which determine the roots 

2az^ - {2hk + a)zi - 2c - 2a ^ = , i = 1, . . . , , (24) 

l^i ~ 

and the equation which gives the energy in terms of the roots Zj's 

TV ^ 

E^2acY,-- (25) 



=1 



Each set of {zi} determine a QES energy E with the corresponding polynomial (f). 

As an example, consider the j — 1/2 case with = 1 and (f){z) — z — Zi. There are two 
solutions. From eq.(24), one sees that the root Zi satisfies 

2azl - {2hk + a)zi - 2c = , (26) 

which gives two solutions 

, {2hk + a) ± J{2hk + Qi)2 + 16ac ^ ^ 

-1 = ■ (27) 

The corresponding energy is 

1 



= 2q;c 

a 



{2hk + a) T ^(26fe + a)2 + 16ac 



(28) 



For the parameters assumed here, the solution with root z^ = —\z^\ < gives the ground 
state, while that with root Zi > gives the first excited state. The superpotential Wi is 
constructed according to eq.(15) 

h'(x) 



Wi{x) = Wo 



h{x) — Zi 
8 



—ax 



a 



+ b + k . 



(29) 



—ae 



1 + |e' 



.ax 



This gives the magnetic field 



Bi = -W[{x) 



—aae 



ax 



(1 + l^f |e°^)2 ' 



(30) 



and the potential that allows these two solvable states is 



V{x) 



= Vi-E- . 



(31) 



With this potential, the ground state and the excited state have energy E — Q and E — 



— E~ — ay{2bk + + 16ac, respectively. Our construction, based on the connection 
between SUSY and QES systems, always sets the energy of the lowest energy state to zero. 

This example should convey the general ideas of our construction. QES potentials and 
magnetic fields for higher degree N are obtained in the same manner. We note that even for 
higher values of N the equation (24) still remains an algebraic equation whose solutions can 
always be found albeit may be numerically. But even then the system remains a QES one. 

We mention here that QES systems belonging to Class IV and VI can be considered in 
a similar manner. 



We shall consider one more class of QES potential, namely, class II of Turbiner's clas- 
sification. The general procedure is the same as that apphed to class I. But unlike class 
I, IV and VI, which are called the first type QES problems, class II, III and V belong to 
the second type. In the first type QES problems, N + 1 eigenstates are solvable for a fixed 
potential with a fixed degree N. For the second type, on the other hand, there are N + 1 
QES potentials difi^ering by the values of parameters and have the same eigenvalue of the 
i-th eigenstate in the i-th potential. For our present problem, each potential corresponds to 
a magnetic field. Below we shall demonstrate this using class II potentials. 




Class II 
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The general form of class II QES potential is 



Vn{x, A) = d2g-4ax ^ 2ade-^"^ + [a^ + 2d {bk - a{N + 1))' 
+ {2abk -aa + X) e^^^^ + hi . 



-2ax 



(32) 



The gauge function is 



9i^) = ^e-'"' + -e-"^ - bkx . 
2a a 



(33) 



As mentioned before, the parameters must be so chosen as to guarantee the normalizability 
of the wave function. For definiteness, we assume a,d > and bk — b + k < 0. 
Letting z = h{x) = exp{—ax), the equation of 4>{z) is 



2dz'^ + 2az + 2bk - a] '^^^ 
J dz 



.2dNz + -- — 
a za 



cPiz) = . (34) 



The differential operator in eq.(34) can also be written as 



N 



Tn = -aJ^J- + 2dJ+ + 2aJ° 
The energy E and the parameter A are given by 

£; = , 



a[- + l]-2bk 



J + constant . 



(35) 



TV 



A = a {2bk -ol)^ — , 



(36) 



i=l 



where the Zj's are to be solved from the Bethe ansatz equations 



dzf + azi + bk 



a 



a 



0, i = l,...,Ar. 



(37) 



The required magnetic field is again given by the roots z^s through eq.(15). For N — 0, one 

has A = 0. 

So far everything appears to be the same as for class I. The main point to note is that 
Vn{x, a) is a function of the parameter A as well as N, and A is determined from the roots 

{zi, Z2j ■ ■ ■ , zn}. For a fixed N there are + 1 possible sets of the roots. Therefore one can 
construct A*" + 1 possible potentials vl^\x) {m — 0,1, . . . , N) for eq.(6) according to 



(m) 



N 



N 
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= Vn{x, A^'")) -E , E^O . (38) 

Here A^"*) is the parameter evaluated using the m-th set of roots of the Bethe ansatz equations 
in eq.(36), and VT^' is obtained from eq.(15) using the same set of roots. We recall here 
that in class I discussed previously, the superpotential Wjy was calculated using the set of 
roots which gives the lowest energy, but here all the + 1 sets of roots have to be used. For 
each potential V^y"^ only one state is solved (by the m-th set of roots) with the same energy 
E = 0. And each potential y^*") corresponds to a magnetic field B^^^ = —{wj^^)'. For the 
family of potentials, only some (generally one) potentials have ground state solved, while for 
others the solvable state is an excited state. In other words, we have a family of magnetic 
fields for which Pauli equation is solvable for one level with the same energy. 

To illustrate these points, let us now give two examples. Consider first the case for A^" = 0, 
giving only the ground state ip{x) = const, x exp{—g{x) with g{x) being given by eq.(33). 
The energy of this state is £^ = 0. In this case the parameter A is A = 0, and the QES 
potential that gives rise to this solvable ground state is 

V{x)^vi'^ ^Vo{x,X^O) , (39) 

which according to eq.(38) is generated by the superpotential 

Wl^'\x) = g'{x) 

= -de-'"^ - ae-°^ -b-k . 

The corresponding magnetic field is 

So = -Wi'\x) 

= -2dQ;e-'°=^ - aae-'*^ . (40) 

Now we come to the case for N — 1. The QES wave function is t/j — {z — Zi) exp{—g), 
and the energy of the state is £^ = 0. The root zi is solved from the Bethe ansatz equation 
(37) 

dzl + azi + bk-^ = . (41) 
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We recall here that we have assumed a,d> and bk — b+k < 0. Eq.(41) gives two solutions 

(0 1) -a^Ja'' + 4d{\bk\ + a/2) 

2d ' ^''^ 

where z^^^ (z^^) corresponds to the solution given by the minus (plus) sign. Since z[^^ < 0, 
the state ijj constructed with zf^^ is the ground state, while that with z^^ is the first excited 
state. According to eq.(38), the ground state is the only QES state for the system with 
potential v/"^ = Vi{x,X^^^) and the excited state is the only QES state for the potential 
Vi^^ — A^^''), where the parameters A*^°'^) are given by 



^(0,1) ^ Q:(2|bfc| +Q;) 



^(0,1) 

^1 



(43) 



These two potentials are generated by the superpotentials 



h{x) - 4"^ 



= - ae-°^ + b-k, m = 0, 1 , (44) 



1 - 4"^e"^ 



with the corresponding magnetic fields being 

B\ ' = —Wl ' \x) 



z\ ' e 



ax 



= -2dQ;e-'°^ - aae-°^ - a' , m = 0, 1 . (45) 

.(0) 



(1 - zr'e^^y 

The point to note is that the energy of the ground state for the potential Vi ''' and that of 
the first excited for v}^^ are both equal to zero, i.e. E — 0. 

The case for class III and V are the same as the present one. We shall not repeat the 
arguments here. 

3 Magnetic field in symmetric gauge 

We now consider the same problem in the symmetric gauge 

A, = yf{r) , Ay = -xf{r) , (46) 
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where r — x -\- y . The magnetic field — B is then given by 



B{x,y)^-2f{r)-rf'{r) . 



(47) 



The Pauh Hamiltonian is 



(48) 



Here is the z-component of the orbital angular momentum, and /' = df /dr. We assume 
the wave functions to have the form 



(49) 



with 



V 



(50) 



i?2(r)e^("*+^)^ 

with integral number m. The function ipmii^, ^) is an eigenfunction of the conserved total 
angular momentum — + (7^/2 with eigenvalue J = m + 1/2. The components Ri and 

i?2 satisfy 

d^ 



dr'' 



1 1 

77? — — 

+ rV' - 2/ (m + 1) - r/' + ^ 



i?i(r) = ER^{r) , 



(51) 



and 



dr'^ 



+ r^f - 2fm + rf + 



, (m+1) 



2 _ 1 
4 



i?2(r) = £;i?2(r) , 



(52) 



where E is the energy. 

The gauge function g{r) for eq.(51), which accounts for the asymptotic behaviors of the 
system, has the general form 



/r 
pf{p)dp- ^Inr , 7=|m| + l/2. 

The corresponding superpotential W{r) is 



(53) 



W{r) = g\r) 
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= r/(r)-J. (54) 

One can check that the potentials in eq.(51) and (52) are given by and respectively, 
where — W'^ T W for positive m > 0. This again shows the SUSY structure of the Pauli 
equation. Hence, the procedure presented in the last section can also be applied in this case 
(for m > 0). 

Our task is to find the form of /(r) such that the Pauli equation is QES. It is seen that 
both eq.(51) and (52) are in the Schrodinger form. So one could try to find /(r) that would 
fit eq.(51) and (52) into the forms classified in [5]. As before, we shall only consider the 
upper component Ri. The lower component can be obtained by SUSY. We found that there 
exist three forms of magnetic fields which make the Pauli equation QES. These three forms 
of magnetic fields give QES potentials that belong to class VII, VIII and IX of Turbiner's 
classification. Below we shall discuss only the case for class VII. The other two classes can 
be considered in a similar manner. 

Class VII 

By inspection one finds that if we let /(r) = /o('^) — cif^ + 6 (a > 0, 6 are constants), 
then eq.(51) belongs to class VII of Turbiner's classification with N — 0. With this form of 
the function /, the potential in eq.(51) is 

V{r) = aV^ + 2a6r^ + [b^ -2a{m + 2)] + 7 (7 - 1) r'^ - 26 (m + 1) . (55) 

The magnetic field is Bq — —Aar"^ — 2b. The general potential in class VII has the form 

VN{r) = o?r^ + 2a6r^ + [6^ - a {AN + 27 + 3)] + 7 (7 - 1) r'^ - 6 (27 + 1) . (56) 

Comparing cqs.(55) and (56) one concludes the potential (55) does belong to class Vll with 
N = and for m > 0. 

As in the asymmetric case, we assume Ri — exp(— g'(r))0, then satisfies the same 
equation (9) with all the derivatives now being with respect to the variable r instead of x. 
With the choice 

— ^^"^ + ^'"^ — 7lnr , (57) 
14 



one can check that V{r) in eq.(55) is generated by Wo(r) = g'{r) in the form V — Wq — 
Wq. Hence the method used in the asymmetric gauge can also be apphed here to generate 
magnetic fields which allow for QES potentials with higher weight. To proceed, we need to 
obtain the Bethe ansatz equations for 0. 

Letting z = h{r) = r^, eq.(9) becomes 



d^chiz) 



\az^ + Abz-2 (27 + 1)1 - [AaNz + E] (l)(z) = . (58) 



dz^ 

In terms of the sl{2) generators J^, J~ and J°, the differential operator in eq.(58) can be 
written as 

Tvn = -4 J° J- + 4aJ+ + 46J° - 2 (A^ + 27 + 1) J" + constant . (59) 

For N — 0, the energy is £^ = 0. For higher N > and 0(r) = IliLii^ — ^i), the function 
/(r) = fNir) is obtained from eqs.(15) and (54): 

M.)=/.„-i|:^. (60) 

For the present case, the roots Zi^s are found from the Bethe Ansatz equations 

2az^ + 2bZi - (27 + 1) - ^ = , i = l,...,N , (61) 



N ^ 



and the energy in terms of the roots Zj's is 

E^2{2^ + 1)J2- ■ (62) 

i=l ^« 

For N — 1 the roots zi are 

^ -b ± ^b^ + 2a(27 + 1) 

= . 63 

^ 2a ^ ^ 

The energies are 

E"^ = 2(6 ± ^62 + 2a(27 + l)) . (64) 

For a > 0, the root z^ — —\zi\ < gives the ground state. With this root, one gets the 
superpotential 

Wi(r) ^ ar^ + br - , _, - - . (65) 

+ Pi I r 
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The QES potential appropriate for the Pauh problem is 



V{x) = - W[ , 
= Vi-E- . 



(66) 



With this choice of the potential, the ground state and the excited state have energy E — 
and E = E^ — E~ = ^^b"^ + 2a(27 + 1). The magnetic field Bi is calculated from eq.(47) 
using the function 



Wr{r) + ^ 

T . 



ar^ + b- 



which gives 



Bi{r) = -4ar^ -26 + 



+ \zi 



4kr 



(67) 



(68) 



(r2+ |^rl)2 ■ 

Just as class I, class VII is also of the first type. On the other hand, class VIII and IX 
belong to the second type. We will not repeat the discussions here. 



4 Summary and Discussions 

In this paper an attempt to give a QES generalization of the result of Aharonov and Casher 
is presented. We have given a general procedure for determining possible (nonuniform) mag- 
netic fields such that the Pauli equation becomes QES based on the sl{2) algebra. This 
procedure makes full use of the close connection between QES systems and SUSY. Of the 
ten classes of s/(2)-based one-dimensional QES systems, we have found that only nine classes 
allow such construction. It would be interesting to extend our procedure to the Dirac equa- 
tion. 

The Pauli equation is supersymmetric owing to the fact that the gyromagnetic ratio is 
two, i.e. g — 2. We would hke to mention that recently it was realized [16, 17] that if one 
changes g to some unphysical values g = 2n (n positive integers), then for magnetic field 
of special exponential and quadratic forms, the generalized Pauli equation could possess a 
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new type of supersymmetry [18, 19, 20]. This is the nonhnear generahzation of the usual 
super symmetry, and is given the name "nonhnear holomorphic supersymmetry" in [16, 17, 
20], or " n-fold supersymmetry" in [19]. It is characterized by a non-hnear superalgebra 
among the supercharges and the Hamiltonian, and the anticommutator of the supercharges 
is a polynomial of the Hamiltonian. The usual SUSY can be viewed as a special case, namely, 
the n — 1 case of the n-fold SUSY. Soon after its discovery, the n-fold SUSY was shown 
to be closely related to quasi-exact solvabihty [19, 20]. For the generahzed Pauh equations 
considered in [16, 17], the weight j = N/2 characterizing the quasi-exact solvabihty of the 
system is given in terms of the number n of the n-fold SUSY and some parameters of the 
system. The authors of [17] found certain duality transformations which mix the number 
n and the parameters to give different values of A^. These duahty transformations thus 
connect different sectors of the generalized Pauli equations. From a mathematical point of 
view, quasi-exact solvability of the generalized Pauli equation is an interesting subject to be 
further explored. It is worth mentioning that the main difference between the generalized 
Pauli equation on a plane considered in [16] and the system considered by us is that in [16] 
the weight j — N/2 is related to the n-fold SUSY hy n — 2j + 1 — N + 1, while in our case n 
is always one (i.e. n = 1) and can be chosen arbitrarily (N = 0,1,2, .. .). Hence when the 
system in [16] is reduced to our case (by setting n = 1), the only QES state that is retained 
is the ground state (corresponding to N — 0). Furthermore, since the number N in our case 
is an arbitrarily chosen number, the kind of duality transformation obtained in [17] does not 
exist in our system. 

Finally we mention a few things about the degeneracy of the energy levels. First we note 
that the Hamiltonians H± are SUSY partners (since they are built from nodeless superpoten- 
tials) and thus shares all the levels of H_ except the zero energy state. Therefore all the 
levels of are doubly degenerate except the zero energy level. This is in agreement with 
the results of ref [1]. We now come to the question of degeneracy of the levels within one 
component, namely, H_. Since in all the cases considered in this paper the magnetic fields 
are nonuniform, so according to [21] the excited states are nondegenerate. This behaviour 
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of the excited states is in contrast to the ground state which is always degenerate with the 
degeneracy depending on the magnetic flux. 
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